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Abstract 

We compute the dynamic structure factors of a dense binary liquid mixture. 
These describe dynamics on molecular length scales, where structural relax- 
ation is important. We find that the presence of a few large particles in a 
dense fluid of small particles slows down the dynamics considerably. We also 
observe a deep narrowing of the spectrum for a disordered mixture composed 
of a nearly equal packing of the two species. In contrast, a few small particles 
diffuse easily in the background of a dense fluid of large particles. We expect 
our results to describe neutron scattering from a dense mixture. 



I. INTRODUCTION 

The behavior of a fluid at large length and time scales is well described by hydrodynamics. 
In recent years hydrodynamics has been extended to molecular length scales [1]. This gen- 
eralized hydrodynamic theory (GHT) was motivated by the desire to describe the neutron- 
scattering spectrum S(k, u) of a dense liquid. S(k, u), which reflects the molecular scale dy- 
namics of the fluid, has a central diffusive peak which dramatically narrows with increasing 
density [2]. GHT provides a physically appealing description of this phenomenon. Struc- 
tural relaxation is very slow on molecular scales in a dense liquid due to close packing of the 
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molecules. Thus density fluctuations decay very slowly via self diffusion at k ~ k where 
S(k) has its maximum. The emergence of this long-lived mode leads to a reduction in the 
half-width of the scattering spectrum. 

In the past few years, much attention has been focused on the study of the dynamical prop- 
erties of dense mixtures [3]. A dense binary mixture is known to be a better glass former 
than a one component liquid. Therefore, an understanding of the dynamics of a dense binary 
mixture is of considerable interest in analyzing the slow relaxation process that ultimately 
leads to the glass transition. 

Recently, the intermediate wave-vector dynamics of a dense binary hard-sphere mixture has 
been studied [4] in terms of the density fluctuations of the two species. The main approxi- 
mations that have gone into this study are the following. First of all, interdiffusion has been 
neglected. Secondly, the wave-vector dependence of the self-diffusion coefficients has been 
neglected. 

A more detailed analysis of a dense binary hard-sphere mixture has been carried out in Ref. 
[5]. Here, the full linear generalized hydrodynamics has been solved and it has been explic- 
itly demonstrated that the density fluctuations are the only slow fluctuations on molecular 
length scales. In addition, both self-diffusion and interdiffusion effects with their wave-vector 
dependence, have been retained. Interdiffusion is a physical process special to a mixture, 
not shared by a one-component fluid. Therefore, inclusion of interdiffusion elicits the special 
characteristics of the dynamics of a mixture. The wave-vector dependence of the transport 
coefficients brings out the nonlocal effects arising from the short-ranged interaction potential 
in a dense liquid. Furthermore, the analysis presented in [5] incorporates certain additional 
static couplings in the extended hydrodynamic equations which had been neglected earlier. 

In the present paper, we use the the density modes obtained in [5] to evaluate the dynamic 
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structure factors of a dense binary mixture. We expect the dynamic structure factors pre- 
sented here to describe the molecular-scale dynamics of a dense liquid mixture measured by 
neutron scattering. 

The paper is organized as follows. In Sec. II we present a theoretical derivation of the 
partial dynamic structure factors of a dense mixture. In Sec. Ill we numerically calculate 
S(k, uj) for a few specific values of the parameters of the mixture and make some predictions. 
Finally, in Sec. IV we end this paper with a few concluding remarks. 



II. THEORETICAL DERIVATION 

The differential scattering cross section of a binary mixture measured in a neutron-scattering 
experiment is proportional to the dynamic structure factor S(k, uj) given by [6] 

S(k, uj) = x 1 b 2 S n (k, uj) + x 2 S 22 (k, uj) 

+2^xTx^bS 12 (k,uj), (1) 

where Xi = ni/{n\ + n 2 ) is the relative number concentration in the mixture of particles 
of component i ■ b is the ratio of the scattering length of the first component to that of 
the second. The partial dynamic structure factors Sij(k,uj)(i, j = 1,2) appearing in (1) are 
given by 

1 f°° 

Sij(k,u;) = — dt e^Fiifat). (2) 

Z7T J-oo 

These are the temporal Fourier transforms of the time correlation functions F^ik, t) defined 
as 

Fy(k,t) = (8nMp)5n ji (t)). 

The brackets () indicate an equilibrium grand canonical ensemble average. Sn^i = 1,2), 
the fluctuation of the number density of the ith component at wave- vector k ^ is given by 
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for a fluid containing Ni of type i particles, a labels the particle number. 



It has been shown in [5] that momentum and temperature fluctuations decay very fast on 
molecular length scales in a dense binary liquid mixture. Therefore the molecular-scale 
hydrodynamic description of such a system involves only the partial density fluctuations 
of the two species, which are long lived on such length scales. These fluctuations are the 
intermediate wave-vector extensions of the long-wavelength heat and diffusion modes. 



It is convenient to consider linear combinations a^(fj, = 1, A) of the number density fluctu- 
ations Sn^ and 5n 2 ^ that form an orthonormal set [5]. We choose 

5n, 
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and 
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■ M ' cos a(k) 
8n. 



a\ Sn ik 
sina(fc)- 



y/niS n (k) 
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^Jn 2 S 22 {k)_ 



Here a(k) is given by 



a(k) = sin 1 



.y/Su(k)S 22 (k)\ ' 

with Sij(k)(i,j = 1,2) representing the partial static structure factors of the mixture [5]. The 
time evolution of the densities and a A ^ is governed by the pseudo-Liouville operator [5] . 
The Laplace transform of the generalized hydrodynamic equations for the number density 
fluctuations of a dense binary liquid mixture may be written in the following compact matrix 
form: 



[zl + A(k 1 ,z)]\a(k,z)) = \a(k,t = 0)), 



(3) 



where \a(k,z)) represents a two component column vector with a^z) and a A ^(z) as its 
entries. a^(z)(i = 1, A) are expressed as 



too 

a rk (z) = I dt e~^\ u 
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(z - L + ) 

for Re(z) > 0, with L + the pseudo-Liouville operator defined as in [5]. The hydrodynamic 
matrix A(k, z) has the following elements: 



A n (k,z) = M u (k,z) - 
A 1A (k, z) = M 1A (k, z) - 



n u (k) ' 
Q u (k)Q Al (k) 



ttu(k) 
= A A1 (k,z), 

and 

A AA (k,z) = M AA (k,z)-^/£. 

Here Mij(k,z) with i,j = 1, A and flij(k) with i = 1, A, I and j = I are given in [5]. The 
solution to (3) can be written as 

\a(k, z)) = R(k, z)[zl + A d {k, ^)]- 1 R T (A;, z)\a(k, 0)). (4) 

In Eq. (4) A d (k, z) = H T (k, z)A(k, z)~R(k, z) is a diagonal matrix with the eigenvalues 



[Au(k, z) + A AA (k, z)\ ± J[A u (k, z) - A AA (k, z)] 2 + 4A lA (k, zf 
Z± = 2 

of the matrix A(k,z) as its diagonal elements. The eigenfunctions \Q a )(a = ±) of 

the hydrodynamic matrix A(k, z) are column vectors with —, a ia(M) an( j 

J J K ' ; y/A 1 ^(k,zy+[A 11 (k,z)-z a ^ 

[zc-Auik^)] ag faeiY entries. These eigenfunctions form the two columns of the 

V^lA(M 2 + [Aii(fc,*)-* a ] 2 

orthogonal matrix H(k,z) which diagonalizes A(k,z). 



Equation (4) can therefore be recast as 

\8n(k,z)) = G(k,z)\5n(k,0)), 
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(5) 



— * 

where \5n(k,z)) is a column vector with 5n^(z)(i = 1,2) as its entries and the matrix 
G(k, z) is given by 

G(k,z) = B-\k,z)R(k,z) 

x [zl + A d (k, ^)]- 1 R T (A;, z)B(k, z). 

Here B(/c, z)\5n(k, z)) = \5a(k,z)). In terms of components Eq. (5) reads 

Srii(k, z) = Gij(k, z)Srij(k, t = 0). (6) 

Multiplying both sides of Eq. (6) by 5rii(—k, t = 0) and taking a thermal average we obtain 
the partial dynamic structure factors 

2 

Su(k, z) = Gij(k, z)Sji(k). 
i=i 

Finally the partial dynamic structure factors appearing in Eq. (2) are given by 

^,,)=h ^ (M = t -^ . 

We insert the explicit expressions for Sij(k,u)(i, j = 1,2) in Eq. (1) to obtain the dynamic 
structure factors S(k, cu), which are in turn proportional to the differential neutron-scattering 
cross sections. 



III. RESULTS 

We present results for a mixture of hard spheres of diameters U\ and o"2(o" 2 > <Ti), number 
densities n\ and n 2 and of total packing fraction rj = | \n\cr\ + n 2 cr|] = 0.46. The diameter 
ratio <7i/(T2 is taken to be 0.7. We vary the concentration x 2 = n 2 /(ni +n 2 ) of larger spheres 
keeping r\ constant. We have chosen the ratio b of the scattering lengths of the two species 
to be 0.7. Finally we consider mixtures of spheres of equal masses. This is done mainly to 
focus on the role of structural parameters in slowing down the dynamics. 
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Figures 1 and 2 show S*(k,u) [i.e., S(k,u) scaled with its value at to = 0] as a function of 
the frequency uo scaled with a time given hy j- = 4 '\A nCT i2£ 1 I w jth n the total number 
density of the fluid, a 12 = 0-1 +°" 2 , ^ 12 the pair correlation function at contact between type 1 
and type 2 spheres and /ii 2 = m^/ (mi +m2) the reduced mass. All plots of S*(k,cu) have 
been made for /c<7i 2 = 27r. We notice that in Fig. 1, the curve for x 2 = 0.01 is significantly 
narrower than that for x 2 = 0.9. This observation can be interpreted as follows. A dense 
liquid mixture consisting of a few large particles suspended in a liquid of small particles slows 
down the dynamics on molecular length scales more easily than a liquid with a few small 
particles in the background of large ones. This is due to the fact that large spheres diffuse 
very slowly in a matrix of small ones. In contrast, the diffusion coefficient of small spheres 
being large, they diffuse faster in a background of large ones. This observation agrees with 
recent molecular dynamic simulations in dense mixtures [7], where the authors find that a 
liquid mixture with a majority of large spheres is easier to crystallize compared to one which 
consists mainly of small spheres. This stems from the difference in the rate of diffusion of 
the two species as mentioned above. 

In Fig. 2 we notice that a compositionally disordered mixture (x 2 = 0.2) characterized 
by nearly equal packing fractions (rji = 0.58r],r] 2 = 0.42^) of the two species, undergoes 
a dramatic narrowing reflecting the emergence of a very slow dynamical process. In con- 
trast, dynamical relaxation is relatively faster in a mixture consisting mainly of large spheres 
(x 2 = 0.9). In a disordered liquid mixture, due to a difference in the sizes of the two species, 
the dynamics is slowed down considerably. This happens because both small as well as large 
spheres can be trapped in a cage formed by the other particles, resulting in the formation 
of peaks in Sij(k)(i,j = 1,2) at well separated values of the wave vector. In contrast, in a 
liquid consisting mainly of large spheres, caging is relatively ineffective due to the ease of 
movement of small particles in the matrix of large ones. 
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IV. CONCLUSION 



We have obtained the dynamic structure factors of dense mixtures on molecular scales by 
confining ourselves to the modes of total density fluctuations and interdiffusion, which govern 
the slow dynamics on such length scales. The relaxation of density fluctuations on molecular 
scales depends more strongly on the structural parameters, such as the size ratio, than on 
the mass ratio of the two species [5,8]. This is a consequence of the crucial dependence of the 
dynamics on the rigid static structure of a dense liquid on molecular scales. Therefore, we 
have considered mixtures of equal masses but different sizes. We have confined ourselves to 
mixtures of not too disparate scattering lengths in order to have comparable contributions 
from the two species to the total S(k,cu). 

Comparison with earlier work: Our present analysis is more complete compared to previous 
work [4] in that it takes into account effects of interdiffusion, wave-vector dependence of 
transport coefficients, and incorporates a few additional static couplings in the extended 
hydrodynamic equations which had been neglected earlier. We have checked that our re- 
sults converge to that of [4] in the appropriate limit. We would like to point out that the 
half-widths evaluated from our theory differ from those predicted by [4]. For instance for 
x 2 = 0.5 we find that the half-width predicted by [4] is about 1.5 times that obtained from 
the present theory. This observation points to the fact that inclusion of the additional effects 
incorporated in the present analysis makes a significant difference to the dynamic structure 
factor of a dense mixture on molecular scales. Whether the present analysis improves over 
the earlier theory can only be settled by future neutron-scattering experiments. 

As mentioned in Sec. Ill, we have demonstrated that our quantitative evaluation of S(k, u) 
confirm the qualitative observations made in [7] regarding the difference in the dynamics 
of slowing down between a mixture consisting of a few large particles and one containing 
a few small particles. Earlier researchers [6] had focused on neutron scattering in dilute to 
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moderate density (77 << 0.46) fluid mixtures. They had mainly concentrated on unusual 
sound propagation in such systems. The regime of density considered here is quite different. 
We have looked at high density mixtures where structural relaxation is important and the 
sound modes are overdamped. Our results are applicable to dense liquid mixtures below the 
glass transition packing fraction. 
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